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ABSTRACT 

By logging encounters between planetesimals and planets we compute the distri- 
bution of encounters in a numerically integrated two planet system that is migrating 
due to interactions with an exterior planetesimal belt. Capture of an irregular satel- 
lite in orbit about a planet through an exchange reaction with a binary planetesimal 
is only likely when the binary planetesimal undergoes a slow and close encounter 
with the planet. In our simulations we find that close and slow encounters between 
planetesimals and a planet primarily occur with the outermost and not innermost 
planet. Taking care to consider where a planet orbit crossing binary planetesimal 
would first be tidally disrupted, we estimate the probability of both tidal disruption 
and irregular satellite capture. We estimate that the probability that the secondary of 
a binary planetesimal is captured and becomes an irregular satellite about a Neptune 
mass outer planet is about 1/100 for binaries with masses and separations similar to 
transneptunian planetesimal binaries. If young exoplanetary debris disks host a binary 
planetesimal population then outwards migrating outer planets should host captured 
irregular satellite populations. We discuss interpretation of emission associated with 
the exoplanet Fomalhaut b in terms of coUisional evolution of a captured irregular 
satellite population that is replenished due to planetary migration. 



1 INTRODUCTION 



The majority of irregular satellites about the gas giant plan- 
ets in our Solar system are in eccentric and high inclination 
orbits suggesting that likely originated in heliocentric or- 
bit and were later captur ed into orbit within a planet's Hill 
sphere (see the review bv lJewitt fc Haghighipouill2007l ). Ir- 
regular satellite populations could ha ve at one time bee n 
a source of heavy dust production l|Bottke et al.l l201(t ). 
If irregular satellite populations exist in distant planetary 
systems then exoplanets could be discovered by detect- 
ing scattered stellar light from their associated dust clouds 
iKennedv fc W vatt 2011). The color of the visible light from 
the detected object, Fomalhaut b, in orbit about the nearby 
star Fomalhaut . is consistent wit h that of reflected light 
from the star ( Kalas et all |2008| ) and so does not origi- 
nate from the planet itself but cou ld arise fr om a plane- 
tary ring system (as discussed by iKalas et al. 2008), dust 
cloud associated with ir regular satellites (as proposed by 
I Kenned V fc WvatO I2OIII ) or circumplanetary disk akin to 
the birth site of the G alilean satellites (e.g., as explored by 
ICanup fc WaTdll2002h . 

We focus on extrasolar planets that are in proximity 
to a debris disk as this disk could host planetesimals that 
could be captured by the planets in the system, becom- 
ing irregular satellites. Neither the Fomalhaut system or 
the HR8799 system are likely to have strongly scattered or 



thick plane tesimal disks, the Fomalhaut system (|Kalas et al.l 
l2005l . 120081 ') because the disk is thin (wi th disk aspect ra- 



tio onlv fe/r 0.013 jOuillen et al.ll20"07l '). and the HR8799 
system (|Marois et al.l 2010l ) because the planets are closely 
spaced and the time u ntil the planet's cross each others 
orbits is probably short iGozdziewski fc Migaszewskil l2009l : 
iFabrvckv fc Murrav-Clay||2010l : iMoore fc QuillenlbOlll l im- 



plying that the planets have not yet experienced close en- 
counters. The Fomalhaut and IIR8799 systems can be con- 
sidered to be analogs to our solar system prior to the late 
heavy bombardment epoch of our solar system that in- 
volved close approaches between the giant planets (the 'Nice' 
model: Tsiganis et al.,.2005 ) and subsequent scattering of the 
planetesimal disk. 

In our Solar system, the irregular satellites were likely 
captured during planet /planet close approaches dur ing the 
late heavy bombardment era l|Nesvornv et al.ll2007l ). How- 
ever, the recent discovery o f a la rge transneptunian bi- 
nary population ||No11 et all I2OO8I ) has motivated study 
of the exchange reaction mechanism for irregular satellite 
capture ([Ag nor & Hamilton 2006; Vokrouhlicky et al. 200^: 
Phil pott et al.ll2010l : iNogueira et al.ll201ll '). Irregular satel- 
lites in orbit around the giant planets prior to the epoch 
of late-heavy bombardment were likely lost durin g close ap- 
proaches between planets (|Nesvornv et alll2007l l. Little is 
known about possible irregular satellite populations that 
could have existed prior to the late heavy bombardment era. 
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Both exoplanets Fomalhaut b and HR8799 b are located 
just interior to dusty disks that are generated fr om collisions 
of planetesimals and so are called debris disks (|Marsh et alj 
l2005l : ISu et"al]|2009l l. Planetesimals originating from these 
disks can be scattered by the outer planet into orbits that 
cross the orbits of interior planets. The interaction between 
planetesimals and planets allows the outer planet to mi- 
grate outwards into the planetesimal disk via the exchange 
of angular momentum between planet s and planetesimals 
l|Fernandez fc id Il984l : llda et all [20001 ) . As the outermost 
planet migrates outwards the population of planet crossing 
pla netesimals can be continua lly replenished. 

[Kennedy fc WvafrS (|201lh showed that a dust produc- 
ing irregular satellite population could be long lived, how- 
ever they did not explore mechanisms for irregular satel- 
lite capture about exoplanets. If the planetesimal disks in 
these systems contain binary planetesimals, then the ex- 
change reaction mechani sm for irregular satellite capture 
l|Agnor fc Hamiltonll2006l 'l is a viable mechanism for the cap- 
ture of irregular satellites about planets from the popula- 
tion of planet-orbit crossing planetesimals. Currently, it is 
estimated that binaries account for 30% of transneptunian 
or Kuiper belt objects (KBOs) with inclinations < 5°, and 
~ 5% of the rest of the KBOs (jNoU et al.ll2008l ). If binary 
planetesimals are common in extrasolar planetesimal disks 
then the migration of an outer planet into a planetesimal 
disk would cause binary planetesimals to cross the orbits 
of giant planets in the system, facilitating the capture of 
irregular satellites. 

A planetesimal binary is disrupted inside the Hill 
sphere of a planet when the distance to the planet is 
sufficiently small that the planet's tidal force overcomes 
the gravitational attraction of the binary. The disrup- 
tion can leave an object in a bound orbit around the 
planet if the velocity of the binary components is suffi- 
ciently large and oriented so as to decrease the energy 
of the lower mass pla netesimal when the binary is dis - 
rupted (|Agnor fc Hami lton 2006: IVokrouhlickv et aklbOOSl ). 
The probability of irregular satellite capture via binary ex- 
changes depends b oth on the statistics of c l ose approaches 
and their energies JVokrouhlickv et al.ll2008l : iPhilpott et al] 
l20ld : iNogueira et aklboTlir In this paper we consider the 
statistics of close approaches between planetesimals and 
planets in systems containing two migrating planets. Our 
goal is to better understand mechanisms for irregular satel- 
lite capture that might operate in young exoplanetary sys- 
tems and prior to the epoch of late heavy bombardment in 
our Solar system. 



2 CLOSE ENCOUNTER DISTRIBUTION FOR 
A MIGRATING TWO PLANET SYSTEM 

We numerically study systems containing two planets that 
are in proximity to a planetesimal disk that is exterior to 
the two planets. We restrict our study to low eccentricity 
and low inclination systems. We first discuss our numerical 
integrations. We compute the distribution of the first close 
encounter that would disrupt a planetesimal binary. We then 
consider the velocity distribution of these encounters. 



2.1 Numerical Integrations 

We ha ve modified the hybrid symplectic integrator QYM- 
SYM (|Moore fc QuillenI Boill ') to record properties of all 
particles that approach within a Hill radius of a planet. 
The integrator runs at two levels, an outer level where all 
particles are integrated with a symplectic integrator with a 
fixed timestep, r. When a particle, or particles pass within 
a planet's Hill radius the encounter is integrated more care- 
fully with an adaptive stepsize N-body integrator. Here we 
use the word encounter to denote a trajectory that passes 
within the Hill radius of a planet, not a collision with a 
planet. For every particle that passes within a planet's Hill 
radius in each larger timestep, r, we use the adaptive step- 
size finer integration to record the minimum distance to the 
planet, and at that location we record the relative energy 
per unit mass 

= — 2 

where v, Vp are the velocities of particle and planet, r, are 
the positions of particle and planet, respectively, and Mp is 
the mass of the planet. For _Ep > we compute a velocity 

Voo = ^/2Ep (2) 

corresponding to the velocity of the particle distant from 
the planet and with respect to the planet, were the particle 
and planet isolated and not orbiting the central star. This 
velocity can be used to estimate the probability that an ex- 
change reaction with a binary pl anetesimal can occur leav - 
ing behind a bound satellite (e.g. IVokrouhlickv et al.ll2008l ) 
though a more accurate calculation uses the Jacobi constant 
(Philpott et al. 201^ o r dire ctly integrates orbits of all four 
bodies (Nogu eira et al.ll201ll ). Here we do not integrate the 
four body problem for a planetesimal binary moving in the 
gravitational field of a planet in orbit about a star. However 
we do integrate the trajectories of single planetesimals in the 
gravitational field a planet in orbit about a star and using 
these trajectories we estimate the probability that a binary 
would disrupt near the planet. 

We work in units of the mass of the central star M, , 
with distance given in units of the initial semi-major axis 
of the innermost planet, ai, and time in units of the initial 
orbital period of the innermost planet. The masses of the 
planets are Mi,M2 and semi-major axes a\,a2. The initial 
semi-major axes of the planets are ai = 1 and = 1.4. The 
initial planet inclinations and eccentricities were set to zero. 
The initial mean anomaly of the two planets were randomly 
chosen. 

The planetesimal disk properties are identical for all 
integrations. The initial planetesimal disk is comprised of 
A'^ = 8192 objects of mass m — Vd"'^ that are distributed 
between semi-major axes of amin ~ 1.6 and Umax = 2.5. 
The distribution of planetesimal semi-major axes is flat with 
probability independent of a within amin and amax- The 
initial eccentricity and inclination distributions were chosen 
using Rayleigh distributions with the mean eccentricity e 
equivalent to twice the mean value of the inclination i and 
i — 0.01. The initial orbital angles (mean anomalies, lon- 
gitudes of pericenter and longitudes of the ascending node) 
were randomly chosen. The total mass of the disk is nearly 
10"'^ and large enough that the outer planet can migrate 
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a substantial distance through the disk in a few thousand 
orbits making it feasible to run a number of similar integra- 
tions. 

Our integrations primarily have inner planet with Mi = 
10"'^ or Ml = 10~*, however a few were done with interme- 
diate or larger masses. The properties of the numerical in- 
tegrations are summarized in Table [1] Each simulation was 
integrated for at least P = 1900 orbits with orbital periods 
are measured in units of the inner planet's initial orbital 
period. The outer planet migrated outwards in all simula- 
tions, however, in some simulations (those with low mass 
outer planets) the outer planet first migrated outwards and 
then reversed direction. As a consequence we consider the 
statistics of encounters only during the portion of the in- 
tegration when the outer planet migrated outwards. The 
length of time used to measure properties of encounters is 
also reported in Table [1] in orbital periods of the inner- 
most planet. Final positions of the planets aij and a2j 
and the total number of orbit crossing planetesimals, Nc, 
are also listed in Table [1] The change in total energy due to 
numerical integration errors at the end of the integrations 
was |A_E/_Eo| < 10~'* where Eq is the initial energy (see 
[m oore fc Quilled [ioTH for a discussion about the accuracy 
of this integrator) . The simulations were run using NVIDIA 
graphics cards from the GT200 architecture (the Quadro 
FX5800, or the GTX 285) or using Fermi Class Tesla C2050 
cards from the GFIOO architecture. All of these are capable 
of computing in double precision. 

Planetesimals originating in the outer disk are first scat- 
tered by the outer planet. Afterwards they can cross the or- 
bits of both planets. The interaction between planetesimals 
and planets allows the outer planet to migrate outwards into 
the planetesimal disk via the exchange of angular momen- 
tum between planets and planetesimals (|Fernandez fc id 
1 19841 ). While planetesimals are scattered by planets they can 
have pericenter near that of a planet's semi-major axis. On a 
plot of semi-major axis versus eccentricity these can be seen 
as scattering surfaces associated with each plan et, as seen in 
prev ious integrations of planet migration (e.g.. [Kirsh et all 
12009 ). To illustrate this. Figure [T] shows the semi-major axis 
and eccentricity distribution at the end of the integration 
with Ml = 10"^ and M2 = 10"*. At this time approxi- 
mately half of the planetesimals have experienced encoun- 
ters with a planet. The other integrations look similar but 
with different fractions of particles in each scattering surface 
and differing numbers of particles in between the planets in 
orbits that cross the orbits of the two planets. 

2.2 Trajectories of planet crossing particles 

Because we have logged all particle/planet encounters dur- 
ing the integrations we can examine trajectories of individ- 
ual particles. In Figure[2]we show an example of the trajec- 
tory of a disk particle from the simulation with Mi — 10~^ 
and M2 = 3 X 10""* that was ejected near the end of the sim- 
ulation. We identify ejected particles based on their energies 
(computed with respect to the central star); those with posi- 
tive energy have been ejected from the system. In this figure 
we show the particle's semi-major axis, a, pericenter, g, and 
apocenter, Q, as a function of time. We also show simulta- 
neously distance between particle and each planet and the 
velocity Voo computed with equation ^ as described above. 



0.8 




0.5 1 1.5 2 2.5 3 3.5 4 



a 

Figure 1. Semi-major axis and eccentricities of particles at the 
end of the integration with planet masses Mi = 10~^ and M2 = 
10"*. Solid circles are the planets. Scattering surfaces are seen for 
both planets at eccentricities setting pericenters and apocenters 
near the planets. The outer planet has migrated more than half 
way through the planetesimal disk. 

Encounters are only logged when a particle passes within a 
Hill radius of a planet so Voo and r are only plotted during 
encounters. 

When the particle first crosses the orbit of the outer 
planet it is in a low eccentricity orbit. The velocity of the 
particle with respect to the planet is not large and so en- 
counters can be gravitationally focused. We see in Figure [2] 
that a slow and close approach to the outer planet happens 
before close approaches to the inner planet. The eccentricity 
of the particle must increase before its orbit can cross those 
of both planets. The velocity difference between particle and 
planet is larger when there is a large difference between par- 
ticle and planet semi-major axis and eccentricity. The inner 
planet is larger than the outer one in this integration and so 
the inner planet exchanges more energy with the particle. 
The trajectory shows a series of distant approaches to the 
inner planet that slowly increase the energy of the particle. 
Ejection is not due to a single close encounter with the in- 
ner planet. If this particle were a binary planetesimal then 
it would likely disrupt during one of the slow and close ap- 
proaches with the outer planet just after it becomes orbit 
crossing rather than later on when its semi-major axis is 
high and the relative velocity between planets and particle 
higher. The trajectory shown in Figure[2]also illustrates that 
the order of close encounters with the planets is important. 
A fast but close encounter could disrupt a binary without 
allowing a capture to take place, hence both the order of 
encounters and the velocities of disrupting encounters must 
be taken into account to predict the probability of satel- 
lite capture. In the next subsection we compute histograms 
of disruption locations and velocities during disruption for 
each planet. 

2.3 Histograms of disruption locations 

A planetesimal binary composed of two masses m\,m2 and 
separation a_B in the vicinity of a planet with mass Mp dis- 
rupts at a tidal disruption radius 
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Table 1. Integrations 
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This table lists properties of each numerical integration. First and second 
columns give the masses of the planets in units of the stellar mass; Mi, M2. 
The time P in orbital periods of the innermost planet (based on its initial 
orbital period) used to compute statistical properties of encounters are listed 
in column 3. The semi-major axes of the planets after P orbital periods are 
listed as a\j,a2j (columns 4,5). The total number of unique planet orbit 
crossing planetesimals identified is listed in column 6. Column 7 shows the 
log of the probability that the inner planet would capture the secondary of 
a planetesimal binary due to a binary exchange reaction for a binary that 
disrupts at a normalized tidal disruption radius R^^ = 0.1 (as discussed in 
section 3). This probability is given for planetesimal binary with a primary 
of radius si = 100 km but can be scaled to other planetesimal binaries using 
factors given in eauation l29l The normalized disruption radius, Rj^ is defined 
in equation |4] Column 8 lists the log of a similar probability but for the outer 
planet. The probabilities of capture have been estimated from the encounter 
distributions in the simulations and using equation 1129 | l. 



1/3 



rtd ^ CLB 



(^^) (3) 

V mi + m2 / 

fe.g.. lAgnor fc HamiltorJl2006l ). A planetesimal binary that 
becomes planet orbit crossing will be disrupted the first 
time it comes within a tidal disruption radius of a planet. 
Previous studies have measured distributions of closest 
approaches in different contexts l|Levison fc Duncanlll997l : 
IVokrouhlickv et al.ll200^ but have not measured the distri- 
bution of disruption radii that depend on both planet mass 
and the order of close approaches for the idealized migrating 
system studied here. 

A binary disrupts if it passes within its tidal disrup- 
tion radius of a planet (equation |3]) . The ratio, rtd/Mp^^, 
only depends on binary properties (its separation and total 
mass). For each encounter that approaches within a radius 
rtd of a planet with mass Mp we compute the normalized 
disruption distance 

1/3 

Rtd = ( ^ I (4) 



We now consider the range of possible values for Rtd- 
The planet's Hill radius is 



Mp 



1/3 



(5) 



The normalized disruption distance, Rtd, differs from 
rtd/ r Hp by a ratio of planet semi-major axes. However since 
we expect the ratio of planet semi-major axes is of order 
unity the maximum value of Rtd is approximately 1. The 
minimum value of Rtd occurs when a binary has separation, 
a B , approximately equivalent to the radius of its more mas- 
sive body. This gives a limit 



Rtd < 10" 



1 g cm- 



,10 AU 



(6) 



^lAnit 



where ai^inu is the initial semi-major axis of the in- 
ner planet. A binary that disrupts at Rtd has Rtd ~ 
( ^^y^^^ ) only dependent on binary properties. 
This is why we will measure encounter distances in units 
of Rtd- 



where p is the density of the planetesimals. The range of 
allowed values for Rtd sets the range covered by the x-axes 
of our subsequent figures. 

For each planet, we count the number of particles that 
pass within a given disruption radius. At each disruption 
distance each orbit crossing particle is only counted once. 
A particle is counted at a disruption distance for the planet 
that it first approaches within this distance. The result is 
two histograms (one for each planet) that are functions of 
Rtd and that can be used to predict the disruption radius 
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Figure 2. Evolution of a planetary orbit crossing particle. Bot- 
tom: We show the log of the semi-major axis (purple crosses), peri- 
center (light blue stars) and apocenter (dark blue plus signs) as a 
function of time for a particle that was eventually ejected in the 
integration with planet masses Mi = 10~^ and M2 = 3 X 10"*. 
Changes in particle semi-major axes and eccentricity are caused 
by close encounters with planets. Middle: The velocity with re- 
spect to each planet during close encounters. Green points show 
velocities for encounters with the outer planet and red points show 
those for the inner planet. We note that early encounters with the 
outermost planet tend to have the lowest relative velocities. Top: 
the distance between the object and planet divided by the planet 
mass to the third power during close encounters. Green points 
show distances to the outer planet and red points show those to 
the inner planet. The closest approach at t ^ 200 occurs early in 
the simulation and with the outer planet. This approach also is 
slow and so could allow the secondary of a binary planetesimal 
to be captured into orbit about the outer planet during binary 
disruption. The order of the encounters can be important as an 
early fast close encounter can disrupt a binary making subsequent 
irregular satellite capture impossible even if the individual plan- 
etesimals later on experience slower encounters. 




0.001 0.01 Rtd 0.1 



Figure 3. The fraction of orbit crossing particles (y-axis) that 
pass within a disruption distance of Rtd (x-axis) as defined by 
equation ^ of each planet for four simulations each contain- 
ing two planets. The properties of the binary (mass and separa- 
tion) determine the radius at which it would disrupt or Rtd- Four 
simulations are shown here and all have inner planet with mass 
Ml = 10-^. The outer planet mass is M2 = 0.5, 1, 3 or 10 X 10""' 
depending upon the simulation. The top four and thin curves 
show fractions of particles that approach the outer planets in sim- 
ulations (from top to bottom) for simulations with M2 = 0.5, 1, 3 
and 10 X 10~*. The bottom four and thicker curves show fractions 
for the inner planet for the same simulations. Simulations are la- 
belled by the mass of the outer planet. At each tidal disruption 
distance, each particle is only counted once and for the planet at 
which it first approaches within the disruption distance. Weakly 
bound binaries are preferentially disrupted by the outer planet 
(high curves on the right hand side of the plot). The difference 
between the fraction disrupted by the outer planet and disrupted 
by the inner planet increases as a function of the ratio of inner to 
outer planet mass. 

to be disrupted by the inner planet. More binaries are dis- 
rupted at each disruption radius when the planets are more 
massive. The difference between the fraction disrupted by 
the outer planet and that disrupted by the inner planet is 
largest when the planet mass ratio is large. We expect that 
an outer planet would have a much richer irregular satellite 
population than an inner but more massive planet. 



for binaries of a given mass and separation and to predict 
the planet that is responsible for the disruption. The re- 
sulting histograms for simulations with A/i — and 
M2 = {0.5,1,3, and 10} x 10"* are shown in Figure H 
In these simulations about half of the simulated planetes- 
imals have passed within the Hill radius of a planet. The 
histograms are normalized by the number of planet orbit 
crossing particles, Nc, with this quantity listed for each sim- 
ulation in Table 

Because each particle is only counted one (at each dis- 
ruption radius) we can determine which planet causes the 
disruption for each type of binary. Wide and low mass bi- 
naries are disrupted at large disruption radii whereas tight 
and massive binaries disrupted at smaller disruption radii. 
In Figure [3] we see that wide binaries are most likely to 
tidally disrupt near the outer planet. For simulations with 
similar planet masses, tightly bound binaries are more likely 



2.4 Velocity Distributions 

During tidal disruption the probability that the secondary of 
the planetesimal binary is left bound to the planet depends 
on th e incoming velocity of the binary (|Ae:nor fc HamiltonI 
|200(|). The secondary is more likely to be captured into or- 
bit if the incoming velocity is low. For 5 different disruption 
radii we plot in Figure[3]velocity distributions of the disrupt- 
ing encounters for each of the two planets for the simulation 
with Ml = 10"^ and M2 = 10"* (which is also displayed 
as one of the pairs of curves in Figure [S]). The histograms 
give the fraction of particles with each Voo value (with Voo 
computed using equation [21). These plots have been normal- 
ized so that at each disruption radius the sum of the inte- 
grated velocity distributions for the two planets is 1. The 
curves have been smoothed and particles with computed 
negative Ep neglected from the histograms. During encoun- 
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Figure 4. Velocity histograms for the simulation with planets 
Ml = 10~^ and M2 = 10"" that is also displayed in Figure [S] 
For 5 diflferent disruption radii (Rtd) we plot the distributions 
of incoming velocities, Voo, for disrupting close encounters for 
each planet. The disruption radii are labelled in the key on the 
upper right. Thick lines correspond to the velocity distributions 
for the inner planet and thin lines for those of the outer planet. 
Velocities are in units of the inner planet's initial circular velocity. 
The top dotted curve is that for the inner planet and Rt^ = 0.3 
and the bottom solid curve is that for the outer planet and Rt^i = 
0.3. Encounters with the outer planet are more likely to have 
lower encounter velocities than the inner planet, consequently the 
outer planet is more likely to capture a secondary during tidal 
disruption of a planetesimal binary. 



ters some particles have computed negative Ep and appear 
to be bound when they are not because the motion of the 
planet with respect to the central star has been neglected in 
the computation of Ep. 

We see in Figure [4] that the velocity distributions dif- 
fer for each planet with the mean of the distributions sig- 
nificantly higher for the inner planet than the outer planet. 
Low velocity encounters are much more likely with the outer 
planet than the inner planet. Tight binaries can be disrupted 
by the inner planet but have sufficiently high velocities that 
the capture of the secondary into orbit around the planet 
is unlikely. The curves associated with the outer planet are 
higher at zero velocity than those for the inner planet imply- 
ing that orbit crossing planetesimal binaries are more likely 
to leave a bound object around the outer planet even if they 
are more likely to be disrupted by the inner planet. In the 
next section we estimate the probability of irregular satellite 
capture based on the low velocity region of these distribu- 
tions for our different numerical integrations. 

As seen from Table [T] for most of our integrations the 
outer planet migrates a substantial way through the plan- 
etesimal disk. By dividing integrations into different time 
intervals and comparing the distributions based on encoun- 
ters in these intervals we have checked that the distributions 
based on the encounters (as shown in Figures [3l [4] and sub- 
sequent figures) are not strongly dependent upon time. At 
the beginning of the integrations the outer planet is near 
the inner edge of the planetesimal disk. We have found that 
as long as this edge is sufficiently close to the outer planet 
to allow migration to take place during the integration, the 



encounter distributions are insensitive to the disk edge lo- 
cation. We attempted to determine if the distributions were 
sensitive to disk mass by increasing or decreasing the mass of 
the planetesimals. However we cannot significantly increase 
the number of disk particles without increasing the length 
of computation time. If we significantly decrease the disk 
or planet masses then planet migration can be extremely 
slow. Each simulation took a few days to run, consequently 
we were effectively limited to a narrow region of parame- 
ter space in both planet and disk masses. We compromised 
by fixing the disk mass and then choosing only simulations 
in which the outer planet migrated a significant distance 
through the disk. 



3 EXCHANGE REACTIONS 

We now consider the properties of a planetesimal binary dis- 
tribution that could lead to capture of an irregular satellite 
population in a migrating planetary system. We consider 
a planetesimal binary composed of two masses mi , m2 and 
separation, as, in the vicinity of a planet. During tidal dis- 
ruption in a planet's tidal gravitational field, the change in 
velocity of the smaller body, m2, due to the orbital speed 
about the binary's center of mass 



AV2 



mi 



G(mi + 7712) _ 



mi 



-vb 



(7) 



mi+m2\ OS mi + m2 

(following I Agnor fc Hamiltonlliooil ) where the binary inter 
nal velocity 



Vb 



G(mi -I- m2) 
as 



(8) 



To enable capture, the difference in energy due to the 
velocity change must exceed the initial energy of the 
object or m2V^/2 where Voo approximates the velocity 
(with respect to the planet) before the planet enters the 
planet's Hill sphere (equation [2]) . This condition is easiest 
to satisfy for the seco ndary rather than the primary (e.g., 
iGould fc Quillenll2003l ) hence that is the situation we focus 
on here. 

We define the Hill radius of the binary at its birth semi- 
major axis atirth, 

_ /mi -I- 7712 ,„^ 
rsS = CLbirth I ) • 



3Af. J 

The ratio of the tidal disruption radius to the planet's Hill 
radius 



rtd 

THp 



as \ ( dbirth 
THB 



,1/3 



(10) 



where ap is the semi-major axis of the planet. A binary 
cannot have a separation larger than its own Hill radius. The 
above equation implies that the disruption location in units 
of a planet's Hill radius primarily depends on the binary 
separation in units of its own Hill radius. As mentioned in 
section 2.3, Rtd differs from rtd/rHp by a factor of planet 
radii. Thus the normalized disruption radius Rtd given as the 
X-axis in Figure [3] is approximately the ratio of the binary's 
Hill radius to that of the planet. Transneptunian binaries 
tend to hav e separations of ord er aB/rnB ~ 0.01 (range 
0.002-0.03) l|Grundv et all 120111 ) and so would disrupt at 
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radius (from the planet) of about 0.01 a giant planet's Hill 
radius or with Rm of order 0.01. 

As the capture probability also depends on the binary's 
internal circular velocity, vb, (see equation [7| it is useful to 
consider the maximum and minimum values for this velocity. 
The binary must be bound so that as < rHB leading to a 
minimum binary internal velocity 



vkb 



> 10" 



1 /mi+rri2\- 
[ KL ) 

( 100km) (t 



Lg cm 



M, 



(11) 



(l + q6)^ 



where vkb is the Keplerian circular velocity at the binary's 
semi-major axis (around the star) and qb = m2/mi is the 
binary mass ratio. The above is equivalent to requiring vb > 
vhb where the Hill velocity vhb ~ ^Ikvhb and Qk is the 
Keplerian angular rotation rate (or mean motion around 
the star) at the binary's semi-major axis. When the two 
planetesimals are nearly touching the binary has a maximum 
internal velocity of 



Vb 
Vkb 



ps\a 



(12) 



< 0.003 



lOAU 



lg cm~3 



Mq 



Sl 



100km 



The range of binary internal velocities is restricted by these 
two limits and is not large. The maximum is larger when 
the binary in the outskirts of a planetary system as is true 
for transneptunian binaries. 



3.1 Estimating the Probability of Irregular 
Satellite Capture 

In the above section we showed that a binary planetesimal 
disrupts at a distance from a planet in Hill radii that is pri- 
marily dependent on the binary separation in units of its 
own Hill radius (equation I lOp . A planetesimal binary could 
be tidally disrupted but both objects will remain unbound 
unless the disruption event removes sufficient energy that 
one of the objects becomes bound to the planet. We now con- 
sider the probability that the secondary (from the binary) 
is captured as a function of the binary's incoming velocity, 
Voo, with respect to the planet. We will derive a limit on Voo 
that allows a capture to take place. Using this limit we will 
use our simulations (in which we measure Voo for close ap- 
proaches) to estimate the probability that the secondary of 
a binary planetesimal can be captured into planetary orbit. 

The energy per unit mass of a binary approaching a 
planet 



1/2 



(13) 



where we have only taken into account the gravitational field 
of the planet and neglected that from the star. At the tidal 
radius, rtd, the energy per unit mass can be written 



GMp 
rtd 



(14) 



where v(rttj) is the velocity of the orbit at disruption. For 
the secondary, m2, of the binary to become bound to the 
planet, the change in velocity during disruption must remove 
sufficient energy that the energy, E, becomes negative. The 
change in energy due to a velocity kick at rtd caused by 
disruption 

(v(nd) + Ava)^ 



AE 



v{rtdf , , (AV2)\ ^^, 

— = v(rtd)-Av2+ — - — (15) 

with |At;2| given by equation ((TJ and AE > required 
for capture. Let the angle between the two vectors "virtd) 
and Av2 be 9^. Using the previous equation and this angle, 
for capture we require that 

\Av2\\v(rtd)\cosei, -\ ~ — • (16) 

Because the incoming binary is not bound to the planet, 
the orbital velocity at the moment of tidal disruption, v{rtd), 
must always exceed the escape velocity at rtd- The escape 
velocity at the tidal disruption radius is 



capei^td) — 



2GMp 

rtd 



vb2^ 



hf^h^y. (17) 

Vmi +1X12/ 



We can consider two limits. When the orbit is nearly 
parabolic, the orbit speed at disruption v(rtd) is approxi- 



mately equal to the escape velocity v{rtd) = y ~r~f'' H the 
orbit is highly hyperbolic then v{rtd) is approximately equal 
to Voo- The gravitational focusing factor 



rtdV^ 



'^escape i^^td) 

vZ, 



mi + m2 



(18) 



determines the regime; when this factor is low {Af < 1) the 
orbit is highly hyperbolic. 

It is useful to compare the escape velocity at the tidal 
disruption radius to the velocity kick A«2. Using equations 
(fTT)) and (O the ratio 



T^escape \ f^td ) 



AV2 



ks-if^h^Yfnu+n^). (19) 

Vmi + 7712 / V 77ii / 

Since we expect nearly equal mass binaries the above ratio 
should exceed 1. As Vescape{rtd) is a lower bound on the orbit 
velocity at disruption (and Veacape{rtd > AV2) we find that 
the first term in equation ((16} should dominate the second 
term. Dropping the second term in equation (|16p we can 
approximate the condition for capture as 



cos^„ > 



\/2 



2v(rtd)Av2 ' 



(20) 



Assuming a random distribution of binary orientations, 
we can integrate the fraction of the spherical angle that sat- 
isfies the above condition (equation I20p . We can consider a 
critical angle, 9c, that is the value of 9v that gives equality 
in equation (|20p . Integrating over solid angle (assuming that 
all binary orientations equally probable) the probability of 
capture, pc, depends on this critical angle, 9c, with 

Pc — (21) 

Using equation H20p the strongly hyperbolic orbit (A / < 
1 and viftd) ~ Voo) gives for capture 

Voo Voo 7771 -I- m2 



cos 9 c 



2Av2 



Vb 



2m\ 



(22) 
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The nearly parabolic orbit {Af > 1 and v{rtd) ~ 

Vescapeind)) givBS 

COS 9c ~ r~TX 

Note that cosSc is at maximum 1 if the right hand side of 
either the previous two expressions exceeds 1. Both of these 
expressions primarily depend on the ratio Voo/vb- Capture 
is only likely if cos 6c is less than 1 and so can only take 
place if Voc/vb is small. Since Voo/vb must be small then 
Af oc vb/Voo must be large (see equation II Sp implying that 
the encounter must be nearly parabolic for capture to take 
place. This means we can use equation (|23p and ignore equa- 
tion (|22[) . Furthermore the probability drops from 0.5 to 
over a very small range of values for Voo and mass ratio 
Mp/{mi + 7712). We can estimate the condition for capture 
based on setting the equation (|23p to 1 and determining how 
Voo depends the properties of the binary at this this divid- 
ing line. Setting cos Oc — 1 and inverting equation (|23p we 
find that capture of the secondary during tidal disruption is 
likely (and has a probability of about 1/2) when 

^<(^5^)^l + ..)-2i (24) 

VB V mi -I- 7712/ 

and is nearly zero otherwise. Here the binary mass ratio 
qi, = 1X12 /mi. While tidal disruption only depends on ra- 
dius (from the planet) in units of Hill radii, we see here that 
the probability of capture also depends on Voo/vb- As more 
compact and massive binaries have larger orbital velocities 
they are more likely to leave the secondary as an irregular 
satellite during a close approach (as previous ly discussed by 
lAgnor fc Hamiltonll2006l : iPhilpott et"al]|2010l ') . Here we have 
a quantitative, though approximate condition on V^/vb al- 
lowing capture which we can use to estimate capture prob- 
abilities from the encounter velocity distributions measured 
in our simulations. 

In our integrations we measure velocities in units of 
the Keplerian rotation velocity (about the star) of the inner 
planet {vki), hence it is useful to write our condition for 
capture (equation [24} in terms of Voo/vki or 

< (^ir^y ( (mi + m2)ai,.„.y 

Rather than describe binaries as a function of mass, mi+m2, 
and separation, as , we can describe a binary in terms of its 
mass and the normalized tidal disruption ratio Rtd (defined 
in equation |4|). The condition for capture in equation (|24|) 
or H25p becomes 

Voo ^ / ai,init \^ / Mp\ s / mi + m2 \i 

VKi ^ \ rtd ^^^'^ ^ ^ 

x{l + qbrhhi (26) 




We have estimated that the probability of capture is 
about 1/2 when the incoming velocity, Voo, is less than the 
expression we give in equation (|26|l . Our numerically mea- 
sured distributions of Voo can now be used to estimate the 
probability that a distribution of planetesimals will leave 
behind some captured secondaries. Because only the lowest 
velocity encounters are relevant we can take the velocity dis- 
tributions measured from the simulated encounters, shown 
in Figure |4l and focus on only the number of encounters 
with incoming velocity near zero. Figure [S] plots the frac- 
tion of encounters N{Rtd)/Nc passing within a disruption 
distance that have Voo within to 0.05 (the value in units 
of the innermost planet's Keplerian velocity) for the same 
simulations shown in Figure [S] 

The fraction of encounters with velocities below that 
given by equation (|26p can be estimated by multiplying our 
numerically measured fraction N{Rtd) /Nc integrated over 
values of Voo £ [0, 0.05] by the factor 

0.02 (Rtdy-^ [ Mp/MA ^ 

Here the factor of 0.02, and R'J-^^My^ are those in equation 
(|27p . We convert this normalized fraction to a probability of 
irregular satellite capture using an additional factor of 1/2 
since the probability of capture is at most 1/2. Thus the 
probability of a planetesimal capture (based on our numeri- 
cally measured low velocity fraction of planetary encounters) 
is 



Pcap{si, Rtd) 




To estimate the capture probability we have used a fixed 
range for Voo and corrected for the width of this range 
(rather than integrating to the limit given in equation I26p 
so that at least a few dozen particles are counted at low Voo 
for all planets in all integrations. As the velocity distribution 
decreases with increasing Voo for the outer planet and vice 
versa for the inner planet (see Figure 3)) we have somewhat 
(but not significantly) underestimated the capture probabil- 
ity for the inner planet and overestimated it for the inner 
planet. 

We show in Figure [6] the probability of capture of the 
secondary by each planet, for a primary planetesimal radius 
si — 100 km and low binary mass ratio, qt, and for the same 
simulations shown in Figure [3] and [5] For other values of 
stellar mass, Af», planetesimal density, p, binary mass ratio, 
qb, and planetesimal primary radius, si, the probability can 
be adjusted using scaling given in equation ((29}. 

A comparison between Figures [5] and [5] shows that the 
probability of capture is reduced from that computed with a 
fixed velocity range for more weakly bound binaries (larger 
-Rtd). This follows as more weakly bound binaries have lower 
internal rotational velocities and so cannot exchange as 
much kinetic energy during disruption as more tightly bound 
binaries. We see from Figure [6] that the probability of cap- 
ture about the outer planet is only weakly dependent on 
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Figure 5. Fraction of integrated particles witli incoming velocity 
Voo between and 0.05 (in units of the inner planet's Keplerian 
velocity) for the integrations shown in Figure |3] as a function of 
different normalized binary disruption radii, Rtd- The curves have 
been normalized by the number of orbit crossing particles in the 
simulation, Nc- Thick lines correspond to the encounters with 
inner planets and thin lines to those with the outer planets. In 
all cases slow encounters are more likely with the outer planet 
and so this planet is more likely than the inner planet to capture 
satellites through binary planetesimal exchange reactions. 

the mass of the outer planet. However, the contrast between 
the capture probability for the outer and inner planet is 
higher when the planet mass ratio is larger. When the ratio 
of the planet masses M2/M1 is low, the capture probabil- 
ity for weakly bound binaries (large Rtd) is higher than for 
M2/M1 ~ 1. This is counter-intuitive as it suggests that 
additional weakly bound binaries are captured by the inner 
planet when the outer planet is more massive. The probabil- 
ity of capture by the outer planet is remarkably insensitive 
to the ratio of the planet masses. However it would be non- 
trivial to predict the shape of these probability curves and 
how they depend on planet masses and ratios. 

Our estimated cap ture probability is consistent 

with that estimated by iNogueira et al.l (|201ll ) for the cap- 
ture of Triton via binary exchange reactions during Nep- 
tune's slow migration outwards following close encounters 
by the giant planets within the context of the "Nice" model. 
But as shown here only the outermost planet would tend 
to capture satellites during this slow period of migration, 
so other scenarios are required to account for the similar- 
ity of the irregular satellite po pulations of the giant planets 
l|jewitt fc HaghighipouilboOTl ) and so the majority of irreg- 
ular satellites in our Solar system. 

If dust associated with a satellite population is detected 
in the vicinity of an outer exoplanet but there is no in- 
ner planet detected, then it would be helpful if the ratio 
of satellite capture probabilies could be used to place limits 
on the planet masses. In Figure [7| using different simula- 
tions we plot the ratio of the capture probability for the 
outer planet and that of the inner planet as a function of 
the outer planet's mass. The capture probabilities are those 
at a disruption radius Rtd = 0.1 and for primary planetes- 
imal radius s\ — 100 km and are listed for the different 
simulations in Table [T] Figure [7| shows that the difference 



Figure 6. The probability that a binary with primary radius 
si = 100 km would leave its secondary captured as a function 
of normalized tidal disruption radius, Rtd, for the integrations 
shown in Figures |3] and [5] Using equation l|29| l the capture prob- 
ably can be adjusted for different binary planetesimal masses, 
mass ratios and densities. The probabilities have been calculated 
from the fraction of integrated planetesimals approaching planets 
at low velocities and using factors given in equation I I28I I. Thick 
lines correspond to probabilities for the inner planets and thin 
lines to those for the outer planets. In all cases the outer planet is 
more likely to capture satellites than the inner planet. The proba- 
bility that the outer planet captures a 100 km radius planetesimal 
from a binary with a tidal disruption radius rtd/fup in the range 
0.01 to 0.1 is about 1/100. 

in capture probabilities is sensitive to the outer planet mass 
with lower mass outer planets preferentially capturing more 
bodies. However, the ratio of capture probabilities appears 
to be fairly insensitive to the inner planet mass. We note that 
the two planet masses were necessarily constrained to be in 
a narrow range because of numerical limitations (we did not 
want to integrate more than a few thousand orbits). More 
diverse and challenging simulations would be needed to de- 
termine if constraints on planet masses could be based on 
captured satellite populations. For the range covered by our 
simulations the ratio is about two orders of magnitude sug- 
gesting that outer planets would preferentially be detected 
from dust associated with an irregular satellite population. 



4 COLLISIONAL EVOLUTION OF A 
POPULATION OF CAPTURED 
IRREGULAR SATELLITES 

Planetesimals captured into planetary orbit contribute to a 
dust producing coUisional cascade if they collide with other 
irregular satellites or other objects passing with in the Hill 
radius of the planet. [Kennedy fc WvalrtI ( 201 ll ) primarily 
considered the evolution of a primordial irregular satellite 
population. However if a planet migrates outwards then the 
irregular satellite population can be replenished by recently 
captured satellites. The proximity of Fomalhaut's dust disk 
to Fomalhaut b suggests that this planet could be migrating 
outwards. If so the planetesimal disk can provide a continu- 
ing source of irregular captured satellites around Fomalhaut 
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Figure 7. The ratio of the sateUite capture probabiUty of the 
outer planet divided by that of the inner planet {y-axis) is shown 
as a function of outer planet mass (x-axis) on a log scale. Points 
for each simulation are taken from Table [T] For most situations 
the ratio of capture probabilities is two orders of magnitude sug- 
gesting that outer planets would preferentially be detected from 
dust associated with an irregular satellite population. 



b. The rate that irregular satellites are captured can be de- 
scribed as a mass rate, Ms, the mass in irregular satellites 
captured per unit time. The total rate is integrated over the 
planetesimal distribution, however because the planetesimal 
mass distribution is likely dominated by the most massive 
objects it is simpler to approximate the rate by only consid- 
ering the massive binary population; 



Ms 



rhofBqty 



(30) 



where fs is the fraction of mass in massive binaries with 
typical separations, aB,typ, mass ratios, qtyp, and typical to- 
tal mass, mtyp. Here Ptyp is the probability that a binary 
with these properties would be disrupted and leave a bound 
satellite around the outer planet. The mass rate rho is the 
mass of planetesimals from the disk that is launched per 
unit time into orbits that cross the planets. This rate is de- 
termined by the migration rate of the outer planet or 



rho = 2-KapTjda. 



'^disk f^p 



(31) 



where ap is the semi-major axis of the outer planet that 
is migrating at a rate dp into a planetesimal disk of sur- 
face density E^. The planetesimal disk has mass Mdisk ~ 
2-KapWdisk^d where Wdisk is the width of the planetesimal 
belt. In the ca s e of a single pla net embedded in a disk 
jlda et al.ll200(]| : iKirsh et al.ll2009t ). when there is adequate 
disk material to maintain migration then the migration rate 
is independent of the planet mass with 



^p, single^ 



47rEap 



.M. 



disk CLp 



Wdisk 



(32) 



where P is the orbital period of the planet. We can describe 
planet migration rates in terms of this fiducial value and 



using a unitless parameter rjmig with dp 



rjmig f^p, single 



For comparison, in our integrations the outer planets have 



migration rates ranging from ri„ 



0.2 - 1.5. 



Combining equations pOI 1311 132|l we find 



MsP 

Mr, 



Md 



10" 



6 / Mdisk /Mp \ I Wdisk /cip 



0.1 



A/p/M. 

X I 3 ) rjmigQtyp 



\ 0.2 

fB \ ( Ptyp 



(33) 



(^)( 



0.01 



where we have chosen the fraction of mass in binaries, 
/s, a ppropriate for th e low inclination transneptunian ob- 
jects ()No11 et al.ll2008l ). Because the capture probability de- 
pends on the binary velocity (that has only a narrow range) 
low mass objects are unlikely to be captured. The above 
mass rate then represents few and rare events occurring 
over a moderate time period . As d iscussed by previous 
studies l Estrada fc Mosgueiral l2006l : IPhilpott etal] l2010l : 
: Bottke et a l. 2010) a system with only a few irregular satel- 
lites could still experience significant colhsional evolution. 

We can consider a steady state where the captured 
massive irregular satellites drive a collisional cascade in the 
vicinity of the planet with a dust production and mass loss 
rate Md equivalent to the rate that mass is gained into the 
planet's Hill radius via capture of satellites; Md = Ms- 

The dust production rate can be estimated from the 
opacity of the small dust particles in the system, Td with 
radius Sd ~ lO^m. Smaller dust particles are removed from 
the system due to radiation pressu re or PR drag (see dis- 
cussion bv lKennedv fc WvattI [20T1I ) . The rate of collisions 
is ~ Td^ and so the dust production rate 



Md ~ PdSdrlrinirc) 



(34) 



where rc the radius of the dust cloud and the dust has den- 
sity pd- The luminosity of the dust cloud in units of the 
stellar luminosity depends on the opacity and size of the 
cloud 



Ld 



4ap 



(35) 



where Qa is the albedo. The luminosity ratio is an observable 
quantity and so relevant for the discussion of emission from 
Fomalhaut b. 

Inside the Hill radius of a planet we can describe radii 
in units of the Hill radius ^ = r/rn- The angular rotation 
rate Q.{r)/Q.K ~ where 0,k is that of the planet in 

orbit about the star. We define = rc/rn as the irregular 
satellite cloud radius in units of the planet's Hill radius. The 
luminosity ratio and dust production rate can be written 



Md 



2,1/2 2/3 2o 
PdSdTdt,c P apilK 

\0.2J 



(36) 
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Ld 



= Mp/M, and 
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Here we have chosen parameters iTd,Qa,ic, Mp) based on 
previous scenarios that account for Fomalh aut b's emission 
jKalas et al.ll2008l : lKennedv fc Wvattll201ll ). 

So as to compare to a planet migration rate we show 
the dust production rate in units of the planet's mass and 
orbital period, P, 

Md „ 2^1/2 2/3 2 1 /oo\ 

^ p_i ~ 2npdSdTd^c M apMp (38) 

\o.2j \ 10-^ J Vo.oiy 

f Pd \ f \ ( gp Y ( M^Y^ 
l^lgcm-3 ) l^lO^my \IQQAU) \Mq ) 

We can compare the rate that mass is captured into plane- 
tary orbit fequationl33[) with the dust production rate (equa- 
tion [38]). The above extremely modest dust production rate 
is below the mass that could be captured in irregular satel- 
lites. A coUisional cascade could be maintained if even a 
small fraction of orbit crossing planetesimals are captured 
into an irregular satellite population as the outer planet mi- 
grates outwards. If the collision rate and opacity becomes 
high then small objects could be re-accreted onto larger ones. 
At a sufficiently high opacity a dense and thin ring system 

could form. 

As discussed bv [Kennedy &: WvattI (|201ll ) there is only 
a small range of parameter space for possible irregular satel- 
lite clouds that could account for the optical emission in 
the vicinity of Fomalhaut b. Here we find that this con- 
straint might relaxed if the cloud is replenished by incoming 
captured satellites. These would have to be fairly massive 
(and so rare) to be captured, however their lifetime depends 
not only the timescale that they collide with other object 
bound objects but also with the time it takes them to col- 
lide with the population of orbit crossing planetesimals (in- 
cluding planetesimals that are not binaries). As the planet 
migrates outwards the irregular satellite population would 
build up until there is a balance between mass acquired by 
the system (in the form of captured objects that are then 
fragmented by collisions with either each other or external 
objects) and mass lost that is either removed from the sys- 
tem by radiation forces, or accreted back onto other objects. 
If the opacity is sufficiently high then debris could coalesce 
into a ring. 



5 DISCUSSION AND SUMMARY 

Using numerical integrations of planets interacting with a 
planetesimal belt we have explored the statistics of close 
encounters of planetesimals with planets during planetary 
migration. We have focused exclusively on the problem of a 
two planet system in proximity to a cold (in terms of velocity 
dispersion) planetesimal belt exterior to both planets that 
induces planetary migration. If the planetesimal belt con- 
tains binary planetesimals then close encounters with plan- 
ets can disrupt these binaries. To predict the location of tidal 
disruption (at what radius and which planet), we have taken 
into account the order of close encounters with planets. We 
find that the velocity distribution (for velocities with respect 
to the planet) of tidally disrupting encounters differs for the 



two planets. Lower velocity encounters are only likely with 
the outer planet. As a consequence we predict that the prob- 
ability that irregular satellites are captured due to binary 
exchange reactions is higher, by about two orders of magni- 
tude, with an outer planet than an inner one. We estimate 
that planetesimal binaries similar to those in the Kuiper 
belt would have a probability of about 1/100 of disrupting 
and leaving a bound object about an outer migrating planet. 
Our numer ical estimate is c onsist ent with that recently es- 
timated by iNogueira et al.l l|201ll ) for the capture of Triton 
into orbit about Neptune. In migrating systems we infer that 
only the outer planet could gain a significant population of 
captured irregular satellites and other scenarios would be 
needed to account for irregular satellite populations around 
inner giant planets. We note that the simulations consid- 
ered here were restricted by numerical considerations to a 
narrow range of planet and disk masses, only considered cold 
and exterior planetesimal disks, only considered two giant 
planet systems and situations lacking planetary encounters. 
Future studies can expand the types of systems numerically 
integrated to predict the probability of satellite capture for 
more diverse systems. 

We estimate that the tidal disruption radius in units 
of the planet's Hill radius is approximately set by the 
binary separation in units of its own Hill radius. Thus 
if optical emission from exoplanets is associated with a 
dust cloud generat e d by irregular satellites, as explored by 
[Kennedy fc WvattI (1201 ll ). then there is a relation between 
the size of the dust cloud, the planet mass and the prop- 
erties of the primordial binary planetesimal population. We 
estimate that capture is allowed only if the incoming binary 
velocity, with respect to the planet is only a few times the 
velocity of the binary orbit prior to disruption. We have used 
this estimate and our numerically measured encounter dis- 
tributions to estimate the probability that the secondary of 
a binary planetesimal could be captured into orbit about a 
planet. To improve upon our derived capture probabilities 
it is likely that the four body problem (binary planetesi- 
mal encounter with a planet in orbit about a star) must be 
integrated. 

iKennedv fc WvattI { 201 iD suggested that irregular 
satellites can be captured if Fomalhaut b crosses into the 
planetesimal belt. Here we explore a related scenario. We 
consider Fomalhaut b not crossing into the planetesimal 
belt, but migrating outwards and capturing planetesimals 
due to binary exchange reactions with binaries that have 
been perturbed into planet orbit crossing trajectories. A 
crude estimate of the mass in irregular satellites captured 
per unit time suggests that the mass captured per unit time 
could be larger than the estimated dust production rate if 
a collisional cascade operates among the irregular satellite 
population. Consequently there could be sufficient mass cap- 
tured as irregular satellites to fuel such a collisional cascade. 
In a migration scenario, we find that the outer planet is 
likely to gain two orders of magnitude more irregular satel- 
lites than inner planets. It is tempting to consider this sce- 
nario as accounting for the detection of Fomalhaut b but not 
other, possibly more massive planets that co uld exist interior 
to it in the system (|Kenworthv et al]|2009l ). Further study 
is needed to determine whether a migration scenario for Fo- 
malhaut b that included irregular satellite capture would be 
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consistent with the ~ 0.1 eccentricity of the dust ring and 
planet and the location and shape of the ring edge. 
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